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Section  1 
INTRODUCTION 


For  some  years  there  has  been  a  fjrowinp:  interest 
in  the  implications  of  mesoscale  oceanographic  features  for 
acoustic  propagation.  One  ma.jor  obstacle  to  assessing  these 
implications  had  been  the  absence  of  a  propagation  model 
capable  of  treating  the  range-dependent  oceanography.  With 
the  application  to  underwater  acoustics  of  the  parabolic 
approximation  to  the  wave  equation  (Tappert  and  Hardin, 
1973;  Tappert,  1974;  Brock,  1978)  this  obstacle  was  overcome 
and  several  Investigators  have  estimated  the  acoustic  field 
in  the  presence  of  such  mesoscale  features  as  eddies. 
Indeed,  the  spring  1978  meeting  of  the  Acoustical  Society  of 
America  devoted  an  entire  session  to  this  subject. 

Much  of  the  work  to  date,  however,  has  left 
unaddressed  several  significant  questions: 

•  What  Is  a  useful,  objective  measure  of  the 
acoustic  importance  of  eddies? 

•  For  the  purpose  of  acoustic  studies,  can  a 
parametric  oceanographic  model  of  eddies  be 
used? 

The  first  question  gets  at  the  matter  of  how  to  characterize 
the  importance  of  eddies  to  ASW  systems.  Individual  point 
comparisons  of  the  acoustic  field  with  and  without  an  eddy 
present  can  generate  arbitrarily  large  differences  which  are 
largely  misleading.  Averaged  propagation  may  indicate 
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-significant  differences  where  they  exist,  but  may  obscure 
interesting  differences  in  the  distribution  of  energy  with 
range  and  depth.  In  the  work  reported  here  an  attempt  is 
made  to  define  a  measure  of  difference  which  may  come  closer 
to  being  useful  in  the  ASW  systems  analyses.  The  second 
question  is  of  importance  primarily  to  physical  oceano¬ 
graphers.  However,  a  satisfactory  analytical  eddy  model 
would  provide  a  very  useful  parameterization  for  acoustic 
studies.  For  this  reason  a  recently-developed  eddy  model 
(Henrick,  1978;  1979)  is  examined  here  in  an  attempt  to 
characterize  measured  eddies  from  some  Pacific  data  col¬ 
lected  by  ONR. 

In  the  remainder  of  this  report  the  eddy  model 
is  described  and  applied  to  measured  data.  Following  this 
discussion,  some  acoustic  results  are  summarized  and  com¬ 
pared.  Two  major  conclusions  with  respect  to  the  questions 
raised  above  are; 

•  The  objective  acoustic  comparison  developed 
here  shows  some  promise  of  usefulness,  but 
requires  more  evaluation  before  it  can  be 
claimed  to  fulfill  the  stated  need. 

•  The  eddy  model  is  not  sufficiently  flexible  to 
allow  fitting  the  measured  data  with  enough 
accuracy  for  the  acoustic  studies.  Some 
suggestions  for  its  extension  are  offered. 

Finally,  with  regard  to  ASW  implications  of  eddies,  the 
following  conclusion  is  reached: 
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Except  for  the  case  in  which  the  source  (target) 
is  at  the  center  of  an  eddy  and  the  receiver  is 
outside  it,  the  presence  of  the  eddy  has  a 
negligible  influence  on  the  acoustic  field.  In 
the  one  case  of  interest,  average  level  differ¬ 
ences  of  10  dB  are  observed  at  low  frequencies 
(20  Hz)  which  diminish  at  high  frequencies  (300 
Hz) . 
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Section  2 

EDDY  MODEL  DESCRIPTION 


2.1.  Basic  Equations 

For  horizontal  nonaccelerated  flow  in  a  friction¬ 
less  ocean,  a  balance  exists  between  the  Coriolis  forces 
and  the  pressure  gradient.  This  relationship  is  expressed 
in  (1).  In  addition,  from  the  hydrodynamic  equations,  we 
may  describe  the  vertical  pressure  gradient  as  a  function  of 
density  and  gravitational  acceleration  (2).  Differentiating 
(1)  with  respect  to  z  (3),  (2)  with  respect  to  r  (4), 
equating  the  two,  and  rearranging  terms,  we  end  up  with 
(5).  Typically,  9o/9z  is  much  larger  than  9p/9r,  but  the 
v/ 0  factor  in  front  of  the  9p/ 9z  term  reduces  the  impor¬ 
tance  of  that  term.  Consider  the  following  situation,  in 
which  g,  gravitational  acceleration,  is  approximately  10^  m 
sec“2,  o~i  and  f,  the  Coriolis  parameter,  is  approximately 
10“^  sec“^.  This  gives  us  a  value  for  g/Pf  on  the  order  of 
magnitude  of  lO^.  For  v,  current  velocity,  approximately 
100  cm  sec”l ,  and  if  the  vertical  gradients  are  a  thousand 
times  larger  than  the  horizontal  gradients,  i.e.,  9p/9z  = 
103  9p/9r,  9v/9z  is  approximately  equal  to  10^  -  10®. 
We  may  then  effectively  ignore  the  second  term  in  (5), 
which  then  becomes  (6). 
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fv  -  1  II  (1) 

where 

f  ■  Coriolis  Parameter  *  2nsin4i 
V  ■  Current  Velocity 
<p  *  Latitude 
p  *  Density 
p  =  Pressure 
r  »  Horizontal  Range 


l£  = 

32 


pg 


(2) 


3r3z 


£  -  (pfv)  =  pftt  +  fv^ 


3z 


’ll 

dz 


3p 

3z 


(3) 


3^  P  3  .  3P  (4) 

araz^sr’^^  ®ar 


pf 


fv  l£ 

32 


3p 

3r 


-1 
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2.2.  Henrick  Eddy  Model 

The  eddy  model  developed  by  R.  F.  Henrick  at  RPI 
(Henrick,  1978:  1979)  is  defined  under  the  following  assump¬ 
tions  : 

(a)  The  eddy  is  radially  symmetric  around  a 

vertical  axis:  i.e.,  the  eddy  is  circular. 

<b)  The  eddy  experiences  slow  translation. 

Thus,  translational  effects  may  be  ignored, 
and  the  eddy  may  be  evaluated  for  a  "snap¬ 
shot"  in  time. 

(c)  The  environmental  effects  of  the  eddy 
vanish  at  a  distance  ro  from  its  center. 

(d)  Environmental  effects  vanish  at  a  depth 
Zq. 

(e)  Salinity  is  invariant  in  this  model,  fixed 
at  35  o/oo. 

The  eddy  model  is  governed  by  the  stream  function 
seen  in  (7).  This  function,  X(r,z),  consists  of  two  inde¬ 
pendent  and  separable  terms,  one  a  function  of  range,  r,  and 


the  other  a  function  of  depth,  z.  The  constant,  k,  indi¬ 
cates  eddy  rotational  direction,  with  k2.0  implying  a  warm 
eddy  circulating  in  a  clockwise  manner,  and  k<0  indicating  a 
cold  eddy  circulating  in  a  counterclockwise  manner.  Equa¬ 
tion  (9)  shows  the  depth-dependent  component  of  the  stream 
function  X.  The  constant  B  is  empirically  estimated  through 
a  process  which  will  be  discussed  in  detail  later. 


X(  r ,  z ) 


(7) 


where 


-1 


(8) 


»  First  zero  of  *  3.83 


m  =  Ji(x)  =  .582 
1  in&x 


+  »>  Eddy  rotational  direction 


+  *>  Anticyclonic  (warm,  clockwise) 


•>  Cyclonic  (cold,  counterclockwise) 


F^^^(z)  » 


1+Bz)”*  |cos  I  ^ ^^n(l+Bz) 


sin  [iY^^>£n(l+Bz)]| 


(9) 
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(10) 


(1)  .  2tt 

^nd+Bz  ) 
o 


This  model  uses  dimensionless  parameters  scaled  In 
the  following  manner.  All  horizontal  ranges  are  scaled  by  a 
characteristic  length  scale,  L.  All  depths  are  scaled  by 
the  depth  of  the  ocean,  D.  The  characteristic  length  scale 
Is  a  function  of  0  and  S,  where  S  Is  a  dimensionless  scaling 
constant.  This  scaling  constant  Is  a  function  of  g,  gravita¬ 
tional  acceleration,  Cq,  surface  sound  velocity,  and  the 
Coriolis  parameter,  2  0slni)).  For  example,  the  dimensionless 
horizontal  range  used  In  the  eddy  model  Is  found  by  scaling 
the  dimensional  range,  denoted  by  the  star(*),  by  the 
characteristic  length  scale. 

L  =  DS  L  =  Characteristic  Length  Scale  (m) 

D  =  Ocean  Depth  (m) 

®  (Dimensionless) 

latitude 

-2 

9.8  m  sec 

surface  sound  velocity 
(m  sec”^) 

earth's  angular  velocity 
7.29  X  10“®  sec"^ 


<t»  * 

g  * 

ar 

O 

n  * 
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This  stream  function  Is  subject  to  several  con¬ 
straints  which  result  from  the  initial  assumptions  (a)  -  (d) 
discussed  earlier.  The  stream  function  X(r,z)  vanishes  at  a 
horizontal  range  ro  from  the  eddy  center.  Analogously,  it 
vanishes  at  a  depth  Zq.  In  addition,  the  r  derivative 
3X/3r  and  the  z  derivative  3X/3z  vanish  at  ro  and  Zq , 
respectively.  These  constraints  result  from  the  vanishing 
of  eddy  effects  at  ro  and  Zo.  3X/3r  and  3X/3z  represent 
eddy  velocity  and  density  effects,  which  also  vanish  at  ro 
and  Zo.  3X/3r  is  actually  current  velocity  (11),  and  the 
depth  derivative  of  this  function  appears  in  equation  (12). 
The  value  of  3X/3z,  seen  in  equation  (13),  again  shows  two 
Independent  and  separable  terms,  one  a  function  of  range  and 
the  other  a  function  of  depth.  The  depth-dependent  term 
contains  the  constant  B.  We  now  have  an  expression  for  3v/3z 
in  terms  of  3X,/3z.  Returning  to  equation  (6),  we  see  an 

alternative  expression  for  3v/3z.  Equating  the  two,  inte¬ 
grating  with  respect  to  r,  and  evaluating  at  the  boundary 

conditions  to  determine  the  value  of  the  integration  con¬ 
stant,  we  are  left  with  equation  (14),  where  o  is  density 
inside  the  eddy  and  P*  is  static  density,  i.e.,  density 
outside  the  eddy.  We  now  have  a  means  for  expressing  the 
perturbation  to  the  density  field  caused  by  the  presence  of 
an  eddy. 


k 


(So) 


-B(l+Bz)"^/^(e+^)sina 


where 


a  *  e<in(l+Bz) 

e  *  (13) 


(14) 


One  means  by  which  this  density  field  may  be 
characterized  is  the  equation  of  state  developed  by  Uamayev, 
for  salinity  constant  at  35^/oo*  This  is  seen  in  equation 
(15).  Substituting  this  expression  into  equation  (14), 
we  get  equation  (16),  where  T  is  the  temperature  inside  the 
eddy  and  T«  is  static  temperature,  the  temperature  outside 
the  eddy.  Rearranging  this  relationship  and  redimension- 
allzing  9x/3z  by  Uog/fCo»  where  Uq  is  surface  current 
speed,  equation  (17)  results.  Substituting  the  expression 
for  3x/3z  from  equation  (13),  we  find  equation  (18),  which 
expresses  the  relationship  between  the  unperturbed  and 
perturbed  temperature,  T.  Equation  (18)  is  quadratic  in  T, 
and  may  be  solved  for  T.  Sound  velocity  within  the  eddy  may 
then  be  computed  using  Wilson's  equation. 
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Pg(2)  »|l  +  10 


1^1  +  10“^|28.152-.0735T-.00469T^|j  +  45  x  10"^z  (15) 


H  ^  |l0‘^  |-.0735(T-T*)  00469 (T^-T*^)|j 


(16) 


2  *2  ^  av 

,00469(T^-T  ^)  +  .0735(T-T'  )  »  - £-2 

c  3  Z 


(17) 


0*0  *  -10^  Up 

#Ti  \  _  _  O  O 


.  00469  (T‘“-T  ■“)  +  .0735(T-T  )  = 


m- 


^  .k  J  1-^J-  J  (e  ) 
(  I  1*0  /  o  o' 


|-B(  l+Bz)“^/^(E+^)sina||  (18) 


2.3  Empirical  Estimation  of  Model  Parameters 


Given  a  static  temperature  distribution,  i.e.,  the 
temperature  distribution  in  the  absence  of  an  eddy,  static 
density  may  be  expressed  as  a  function  of  temperature 
and  pressure  (depth),  for  salinity  constant  at  350/oo, 
using  Mamayev's  equation  of  state,  seen  earlier  in  equation 
(15).  We  may  then  generate  a  distribution  of  static  density 
as  a  function  of  depth,  o^Cz),  from  the  temperature  data. 
Henrick  fits  this  density  distribution  to  an  algebraic 
function  for  static  density  shown  in  equation  (19).  The 
Levenberg-Marquardt  algorithm,  a  non-linear  least  squares 
regression  algorithm,  is  used  to  fit  this  functional 
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form  for  static  density  to  the  observed  density  (temper¬ 
ature)  values.  Estimates  of  Pq,  A,  and  B,  which  minimize 
the  residuals,  are  returned.  This  process  may  be  extended 
to  fit  this  functional  form  to  temperature  and  sound  veloc¬ 
ity  observations  directly. 


where 


(19) 


DATA 


The  data  used  in  this  analysis  were  collected 
during  May,  1976  as  a  part  of  OMR's  sh  ips-of-opportun  i  ty 
program.  Six  parallel  ships  centered  at  37^M  traversed  a 
2500  km  swath  from  ISS^E  to  1750W.  The  lateral  spacine 
between  ships  was  approximately  45  km.  The  ships  made 
hourly  XBT  casts,  alternating  between  T-4  and  T-7  XBTs .  This 
resulted  in  a  sampling  density  of  27  km.  The  resultant  temper¬ 
ature  data  have  given  us  one  of  the  best  quasi-synoptic 
pictures  of  mesoscale  variability  to  date.  Texas  AfeM 
University  merged  the  thermal  data  with  deep  ocean  clima¬ 
tology,  to  yield  a  three-dimensional  characterization  of 
this  swath,  in  terms  of  temperature(T) ,  salinity(S),  sound 
velocity(c),  and  density  (^t).  A  schematic  of  the  track  may 
be  seen  in  Figure  1.  Appendix  A  contains  isothermal  and 
isovelocity  sections  for  each  of  the  six  ship  tracks.  It  vs 
evident  from  these  temperature  and  sound  velocity  sections 
that  a  significant  amount  of  mesoscale  variability  is 
present  in  this  region  of  the  Pacific. 

We  have  modeled  most  of  the  eddies  on  these  tracks 
using  the  eddy  model  developed  by  Henrick.  The  eddies  have 
been  characterized  in  terms  of  their  horizontal  (rg)  and 
vertical  (Zq)  extent  of  influence,  as  well  as  the  value  of 
the  parameters  Pq,  A,  and  B.  Table  1  summarizes  the  results. 
Eddies  from  two  of  the  tracks  have  been  selected  as  an 
illustration  of  the  modeling  process.  A  cold-core  eddy, 
located  at  450  km  along  the  track  of  the  Schenectady  ,  the 
northernmost  ship,  and  a  warm-core  eddy  from  the  start  of 
the  track  of  the  Bristol  County,  the  southernmost  ship,  will 
be  used  for  illustrative  purposes. 
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Section  4 
MODEL  RESULTS 


The  cold  eddy  on  the  Schenectady  track  has  an 
effective  radius,  r^,  of  188  km  and  extends  to  a  dt,  th  of 
influence,  Zq,  of  2000  m.  The  backzround  sound  velocity 
conditions  have  been  modeled  using  the  Henrick  model  and  the 
resultant  values  of  Pq,  A,  and  B  may  be  seen  in  Table  2. 
Using  these  parameters,  the  background  density,  temperature, 
and  sound  velocity  profiles  may  then  be  modeled.  Figure  2 
shows  a  comparison  between  the  actual  density  profile  and 
the  modeled  density  profile,  determined  by  the  parameters 
pQ ,  A,  and  B.  The  agreement  is  apparently  quite  good  • 
Figure  3  shows  a  similar  comparison  for  the  data  and  modeled 
temperature  profiles.  These  profiles  also  agree  fairly 
well,  although  there  are  some  differences  in  the  upper  fiOO 
m.  These  differences  are  relatively  small,  however,  with 
the  modeled  profile  differing  from  the  data  on  average  by 
only  a  few  tenths  of  a  degree.  The  modeled  sound  velocity 
profile  is  compared  to  the  data  profile  in  Figure  4.  In 
this  instance,  there  are  some  fairly  significant  deviations, 
once  again  occurring  in  the  upper  600  m.  Note  also  that 
while  the  deep  gradients  on  these  profiles  are  the  same,  the 
modeled  profile  shows  an  absolute  sound  velocity  at  these 
depths  which  is  consistently  greater  than  the  sound  velocity 
seen  in  the  data  profile  at  the  same  depths. 

Additional  information  may  be  seen  in  a  comparison 
of  the  actual  and  modeled  eddy  fields.  Figure  5  shows  the 
modeled  isothermal  structure  for  this  cold-core  eddy. 
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Eddy  Model  Parameters 


Figure  5  Eddy  Model  leothenra  for  Cold-Core  Eddy 


-70  80 

RfiNGE  (ICM) 


Figure  6  shows  a  superposition  of  the  model  isotherms  with 
the  actual  isotherms  seen  along  the  track.  Some  significant 
deviations  may  be  seen.  The  most  significant  deviation 
results  from  the  absence  of  eddy  surface  expression  in  the 
model,  whereas  surface  expression  is  clearly  present  in 
the  data.  In  addition,  the  model  does  not  reflect  the 
asymmetry  seen  within  the  data.  As  one  of  the  underlying 
assumptions  for  the  model  is  that  the  eddy  is  radially 
symmetric,  this  represents  a  fundamental  and  irrecoverable 
difference.  These  differences  are  perhaps  not  as  large  as 
they  appear,  as  the  contouring  process  Itself  is  quite 
sensitive  to  small  changes  in  temperature.  Deviations  of 
several  tenths  of  a  degree,  typically  not  considered  exces¬ 
sive,  will  result  in  major  displacements  of  the  isotherm 
depths  in  regions  where  the  temperature  is  nearly  constant. 
Figure  7  shows  model  isovelocity  contours  for  the  cold-core 
eddy.  Figure  8  shows  the  superposition  of  these  model 
contours  with  isovelocity  contours  from  the  Schenectady 
track.  The  same  differences  between  the  model  and  the 
data  which  were  evident  in  the  isothermal  sections  are  also 
evident  in  the  isovelocity  contours.  Once  again,  we  notice 
the  relative  lack  of  surface  expression  in  the  model, 
whereas  surface  expression  is  clearly  indicated  by  the  data. 
In  addition,  the  model  does  not  account  for  the  random 
effects  which  are  characteristic  of  a  particular  isotherm. 

Turning  to  the  warm-core  eddy  from  the  Bristol 
County  track,  the  effective  radius  for  that  eddy  is  173  km, 
slightly  smaller  than  the  radius  for  the  cold-core  eddy. 
The  depth  of  influence  is  slightly  greater,  extending  to 
2800  m.  Background  sound  velocity  data  have  been  used  to 
estimate  the  values  of  the  parameters  Pq,  A,  and  B,  which 
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may  be  seen  in  Table  2.  The  background  density,  temperature, 
and  sound  velocity  conditions  have  been  modeled.  Figure  9 
shows  a  comparison  of  the  density  profile  generated  by  the 
model  to  the  data  density  profile  for  this  warm-core 
eddy.  As  in  the  cold-core  eddy  case,  the  agreement  between 
the  density  profiles  is  quite  good.  Figure  10  shows  a 
similar  comparison  for  the  background  temperature  profiles. 
Again,  the  deviations  in  the  temperature  profiles  are  more 
apparent  than  those  in  the  density  profiles,  but  the  devia¬ 
tions  averaged  over  depth  are  generally  less  than  several 
tenths  of  a  degree.  The  sound  velocity  comparison  seen  in 
Figure  11  does  show  some  major  differences.  Once  again,  the 
area  of  largest  deviation  is  in  the  upper  POO  meters  of  the 
water  column.  Further,  the  deep  gradient  on  the  the  model 
profile  is  significantly  smaller  than  that  seen  in  the  data 
profile.  This  has  some  serious  implications  for  the  acous¬ 
tical  analysis  which  will  be  discussed  in  detail  later. 
Figure  12  shows  the  eddy  model  isothermal  and  isovelocity 
sections  for  this  warm-core  eddy.  Figure  13,  which  is  a 
superposition  of  these  model  isotherms  and  isovelocity 
sections  with  the  data  isothermal  and  isovelocity  sections, 
again  illustrates  the  lack  of  surface  expression  in  the 
model,  the  inability  of  the  model  to  account  for  asymmetry 
which  is  present  in  the  data,  and  also  the  failure  of  the 
model  to  account  for  random  effects  along  a  given  isothermal 
or  isovelocity  contour. 
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Section  5 

EDDY  MODEL  COMPARATIVE  STATISTICS 


In  order  to  assess  more  accurately  the  dif¬ 
ferences  between  the  modeled  density,  temperature,  and  sound 
velocity  profiles  and  the  corresponding  data  profiles, 
direct  comparison  of  the  difference  in  sound  velocity  at 
discrete  points  along  the  track  is  quite  helpful.  That  is, 
for  every  location  in  which  a  temperature  or  sound  velocity 
observation  is  available,  the  observed  sound  velocity 
profile  may  be  compared,  on  a  point-by-point  basis,  to  the 
modeled  sound  velocity.  Such  a  comparison  for  the  cold  eddy 
is  shown  in  Table  3.  The  mean  sound  velocity  difference, 
seen  at  the  bottom  of  the  table,  is  actually  quite  small, 
less  than  1.2  m  sec“l .  Similarly,  the  mean  temperature 
difference  over  all  depths  is  small,  less  than  .4^C.  This 
is  apparently  quite  good  agreement.  Looking  in  more  detail 
at  the  differences  in  sound  velocity  at  selected  depths,  it 
may  be  seen  that  the  largest  differences  are  indeed  concen¬ 
trated  in  the  upper  600  m  of  the  profile.  This  rein¬ 
forces  the  information  seen  in  the  graphic  display  of  sound 
velocity  profiles. 

Similar  statistics  for  the  warm-core  eddy  may  be 
seen  in  Table  4,  which  once  again  shows  a  mean  sound  veloc¬ 
ity  difference  and  a  mean  temperature  difference  which  are 
quite  small.  The  major  differences  remain  concentrated  in 

the  upper  400  to  600  m,  with  fairly  good  agreement 

^  ■ 

existing  below  those  depths. 


5-1 


SOUND  VELOCITY  DIFFERENCE  (DATA-MODEL)  m  sec-1 
RADIAL  DISTANCE  FROM  EDDY  CENTER 


Depth(m) 

188  km 

155  km 

77  km 

0 

-2.50 

1 

• 

o 

-12.70 

200 

8.04 

7.36 

-  1.19 

400 

1 

• 

o 

o 

• 

o 

1 

-  .49 

600 

-3.16 

-2.70 

-  .02 

800 

-2.29 

-1 .96 

.52 

1000 

-1.52 

-1.29 

.41 

1200 

-0.79 

-  .64 

.46 

1400 

-0.27 

-  .18 

.48 

1600 

.07 

.13 

.47 

1800 

.35 

.38 

.47 

2000 

.52 

.52 

.44 

^Asv 

-  .24 

.05 

-  1.01 

^at 

.25 

.33 

-  .01 

Table  3 

Comparative  Statistics  Cold  Eddy 


SOUND  VELOCITY  DIFFERENCE  (DATA-MODEL)  m  sec'l 


RADIAL  DISTANCE  FROM  EDDY  CENTER 


Depth(m) 

95  km 

65  km 

33  km 

0 

5.18 

1.93 

1.64 

200 

15.79 

17.00 

14.49 

400 

12.03 

12.69 

16.47 

600 

.56 

-  .99 

4.95 

800 

-  .56 

-1.79 

-  .62 

1000 

.15 

-  .81 

-  .10 

1200 

.76 

.03 

.51 

1400 

1.15 

.59 

.92 

1600 

1.36 

.95 

1.15 

1800 

1.44 

1.15 

1.29 

2000 

1.43 

1.23 

1.33 

2200 

1.31 

1.18 

1.25 

2400 

1.10 

1.02 

1.11 

2600 

.79 

.76 

.83 

2800 

.49 

.49 

.56 

3000 

.07 

.07 

.12 

^Asv 

1.77 

1.41 

1.91 

^AT 

.75 

.66 

.79 

Table  4 

Comparative  Statistics  Warm  Eddy 
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Section  fi 
FDDY  MODEL  SUMMARY 


The  eddy  model  developed  by  Henrick  has  been 
examined  in  some  detail.  This  model  provides  an  empirical 
fit  to  the  static  or  backftround  density,  temperature,  and 
sound  velocity  profiles  for  a  particular  ocean  region. 
These  profiles,  and  thus  the  region,  are  characterized  in 
terms  of  the  parameters  A,  and  6.  The  model  then  gener¬ 
ates  ihe  temperature  and  sound  velocity  fields  within  the 
eddy,  giving  us  an  idea  of  the  perturbations  to  these  fields 
which  are  caused  by  the  presence  of  the  eddy. 

An  evaluation  of  the  model  reveals  that  the  mean 
temperature  and  sound  velocity  differences  are  small. 
Statistically,  the  modeled  profiles  are  actually  a  good  fit 
to  the  data.  It  has  been  shown  that  the  largest  differences 
are  concentrated  in  the  upper  600  meters.  These  differences 
are  attributable  to  the  fact  that  the  model  does  not  permit 
surface  expressions.  Further,  the  model  fits  the  background 
conditions  to  a  smooth,  highly  regular  profile.  Thus,  any 
data  profiles  which  are  not  smooth  and  regular  will  differ 
significantly  from  the  modeled  profiles.  These  differences 
appear  primarily  in  the  upper  600  m,  typically  the  region  of 
the  water  column  experiencing  the  most  variability.  This 
poses  somewhat  of  a  problem  for  the  acoustician,  as  most  ASW 
cases  of  interest  have  the  source  located  within  that  region 
of  the  water  column.  The  use  of  this  model  for  acoustic 
analysis,  therefore,  will  result  in  significant  differences 
between  the  observed  and  modeled  conditions.  On  the  basis 
of  this  analysis,  it  appears  as  if  a  model  upgrade  is  in 
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order  before  it  can  be  used  successfully  for  ASW  prediction. 
One  suggestion  for  upgrading  the  capabilities  of  this  model 
might  include  the  use  of  cubic  splines  or  other,  more 
flexible,  functional  forms,  which  will  give  better  local 
approximations  to  these  background  profiles:  also  required 
is  an  extension  to  permit  eddy  surface  expressions. 

Originally,  it  was  intended  that  this  model  be 
used  to  fit  the  data  from  the  six  ship  tracks  and  to 
parameterize  the  eddies  on  these  tracks  in  terms  of  the  the 
constants  Pq,  A,  and  B.  The  modeled  eddies  were  then  to  be 
used  for  acoustic  analysis.  The  differences  between  the 
modeled  and  observed  sound  velocity  profiles  resulted  in 
significant  differences  in  the  acoustic  fields  generated  by 
these  profiles.  Preliminary  acoustic  analysis  of  the  model 
results  showed  significant  deviations  in  the  manner  in  which 
sound  propagates  under  the  modeled  conditions.  For  example, 
the  differences  in  the  deep  gradient  seen  in  the  background 
sound  velocity  fit  for  the  warm-core  eddy  resulted  in  major 
changes  in  the  cycle  distance  of  rays  propagating  under 
these  conditions.  It  was  therefore  judged  that  the  model  in 
its  present  form  was  unacceptable  for  acoustic  analysis,  and 
that  the  most  accurate  assessment  of  the  acoustic  properties 
of  these  eddies  would  result  from  direct  use  of  the  measured 
oceanographic  data.  The  following  acoustic  analysis  there¬ 
fore  is  restricted  to  analysis  of  these  data. 
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Section  7 
ACOUSTIC  ANALYSIS 


7.1.  Introduction 

The  underlying  purpose  behind  the  following 
acoustic  analysis  is  to  assess  the  acoustic  impact  of 
mesoscale  fluctuations  as  a  function  of  varying  source- 
receiver  geometry  and  source  frequency.  To  that  end, 
two  tracks  have  been  constructed  for  use  in  propagation 
modeling.  The  tracks  contain  the  static  or  background  sound 
velocity  profiles  at  their  start.  The  background  profiles 
are  extended  for  185  km.  Immediately  following  this  range- 
independent  185  km  is  a  section  of  track  containing  an 
eddy,  either  the  warm-  or  cold-core  eddy  from  the  Schenec¬ 
tady  or  Bristol  County  track,  respectively.  The  diameter  of 
both  eddies  was  approximately  370  km.  Following  the 
eddy,  there  is  another  185  km  of  background  sound  velocity 
conditions.  We  therefore  have  constructed  two  acoustic 
tracks,  consisting  of  185  km  of  range-independent  conditions, 
followed  by  370  km  of  track  containing  an  eddy  (warm  or 
cold),  followed  by  an  additional  185  km  of  range-independent 
track.  For  comparative  purposes,  a  totally  range-indepen¬ 
dent  track  has  been  constructed,  where  the  background  sound 
velocity  profile  which  exists  at  0  km  is  extended  for  the 
length  of  the  track,  in  this  instance  740  km. 

To  evaluate  the  impact  of  varying  source-receiver 
geometry,  we  have  selected  four  different  situations.  The 
first  has  the  source  located  at  0  km  at  the  start  of  the 
totally  range-independent  track.  The  second  geometry  has 
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the  source  located  at  0  km.  along  the  range-dependent  track; 
l.e.,  the  track  containing  either  the  warm  or  the  cold  eddy. 
The  third  geometry  has  the  source  one-half  of  a  convergence 
zone  closer  to  the  eddy,  placing  it  at  30  km  along  the 
range-dependent  track.  Finally,  the  effect  of  the  source 
located  at  the  eddy  center  along  the  range-dependent  track 
has  been  considered.  The  depth  of  this  source  was  held 
constant  at  30  m  throughout  this  analysis.  The  receiver 
depths  considered  were  95,  450,  700,  and  900  m.  Three 
representative  frequencies  were  selected  for  analysis. 
Those  frequencies  were  20,  300,  and  1000  Hz.  The  Parabolic 
Equation  (PE)  model  has  been  used  to  model  propagation  loss 
(Tappert  and  Hardin,  1973;  Tappert,  1974;  Brock,  1978). 

7.2  Cold  Eddy  20  Hz. 

Transmission  loss  calculations  as  a  function  of 
horizontal  range  (in  km)  for  a  20  Hz  source  located  at 
30  m  as  seen  by  a  receiver  located  at  95  m  may  be  seen 
in  Figure  14.  Transmission  loss  has  been  calculated  for 
each  of  the  source-receiver  geometries  discussed  previously, 
viz.,  a  totally  range-independent  environment,  a  range- 
dependent  environment  with  the  source  located  at  0  km, 
a  range-dependent  environment  with  the  source  located  at 
30  km,  and  for  a  range-dependent  environment  with  the 
source  located  at  the  center  of  the  eddy.  For  clarity  of 
presentation,  the  transmission  loss  curves  have  been  aver¬ 
aged  over  a  horizontal  range  of  one  convergence  zone,  60  km. 
Thus,  only  variations  occurring  on  a  scale  greater  than  60 
km  are  evident  in  these  plots.  In  general,  all  four  source- 
receiver  geometries  considered  here  experienced  similar 
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Figure  14 


levels  of  transmission  loss.  The  track  with  the  source 
located  at  the  eddy  center  experiences  1  -  2  dB  less  loss 
than  the  other  source  receiver  eeometrles.  This  Is  particu¬ 
larly  evident  In  the  first  280  km  of  the  track.  Also  of 
note  Is  the  close  correspondence  between  the  range-dependent 
environment  with  the  source  at  0  nm.  and  the  one  In  which 
the  source  Is  located  one  half  of  a  convergence  zone  further 
Into  the  track,  at  30  km. 

Unaveraged  transmission  loss,  corresponding 
to  the  totally  range  Independent  environment  and  the  range 
dependent  environment  with  the  source  at  the  eddy  center 
from  Figure  14,  may  be  seen  in  Figure  15.  The  conver¬ 
gence  zones  In  the  range-independent  environment  remain 
narrow,  well-defined,  and  clearly  discernable  for  the  length 
of  the  track.  In  contrast,  the  convergence  zones  in  the 
range-dependent  environment  with  the  source  at  the  eddy 
center  broaden  rapidly  and  show  evidence  of  increasing 
multipath  at  the  farther  ranges.  While  the  peak  energy 
levels  in  the  range-independent  case  are  higher  than  those 
in  the  range-dependent  environment,  the  broader  convergence 
zones  in  the  range-dependent  case  result  in  this  environment 
containing  more  energy,  on  average,  than  the  range-indepen¬ 
dent  case.  This  was  shown  in  Figure  14. 

Moving  the  receiver  slightly  deeper  in  the  water 
column,  to  450  m,  produces  transmission  loss  seen  in  Figure 
16.  Here  we  see  that  the  range  dependent  environment  with 
the  source  at  the  eddy  center  shows  clearly  less  loss 
than  the  other  environments  considered.  The  totally  range 
independent  environment  exhibits  the  highest  loss  levels,  at 
times  as  much  as  7  dB  more  than  the  environment  with  the 
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source  at  the  eddy  center.  This  is  also  evident  in  Figure 
17,  which  shows  unaveraged  transmission  loss  for  the  two 
extreme  cases,  the  range->independent  environment  and  the 
range-dependent  environment  with  the  source  at  the  eddy 
center.  Convergence  zones  in  the  range  dependent  case  shift 
and  overlap  early  in  the  track.  Further,  there  are  many 
additional  arrivals  seen  in  the  range-dependent  track  that 
are  absent  in  the  range-independent  environment.  These 
conditions  result  in  the  overall  lower  level  of  propagation 
loss  seen  in  the  range-dependent  environment. 

Figure  18,  which  demonstrates  the  effect  of  a 
receiver  positioned  at  700  m,  near  the  sound  channel  axis, 
shows  similar  findings.  Once  again,  the  range-dependent 
environment  with  the  source  at  the  eddy  center  clearly 
exhibits  less  loss.  A  receiver  located  at  900  m  (Figure  19) 
shows  the  same  trends. 


7.3  Warm  Eddy  20  Hz. 

Corresponding  transmission  loss  figures  for 
a  warm  eddy  may  be  seen  in  Figures  20  -  23.  Figure  20, 
which  shows  the  effect  of  a  receiver  positioned  at  95  m 
in  the  water  column,  shows  a  relationship  between  source- 
receiver  geometry  which  is  the  inverse  of  that  seen  in  the 
cold-core  eddy.  In  this  instance,  the  range-dependent 
environment  with  the  source  at  the  eddy  center  shows  consis¬ 
tently  higher  levels  of  loss  than  do  the  other  environments. 
The  range-independent  environment  in  this  instance  shows 
consistently  less  loss.  These  differences  are  even  more 
apparent  when  considering  a  receiver  located  at  450  meters 
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in  the  water  column  (FiRure  21).  In  this  instance,  the 
difference  in  transmission  loss  between  the  range-dependent 
environment  with  the  source  at  the  eddy  center  relative  to 
the  other  geometries  considered  may  be  as  much  as  10  dB, 
with  the  source  at  the  eddy  center  experiencing  higher 
levels  of  loss  within  the  warm  eddy.  Moving  deeper  in  the 
water  column,  a  receiver  located  at  700  m  (Figure  22)  shows 
the  same  general  trends,  although  the  differences  between 
the  various  geometries  are  not  as  marked.  Figure  23  shows 
transmission  loss  through  a  warm-core  eddy  for  a  receiver 
depth  of  900  m.  Again,  the  relative  levels  of  transmission 
loss  are  consistent  with  previous  observation.  The  environ¬ 
ment  with  the  source  positioned  at  the  eddy  center  shows  as 
much  as  5  dB  more  loss  relative  to  the  other  environments. 

Comparison  of  Figures  20  -  23  shows  the  effect  of 
increasing  receiver  depth  on  transmission  loss.  It  is 
apparent  from  a  comparison  of  these  figures  how  the  distri¬ 
bution  of  energy  changes  as  a  function  of  receiver  depth. 
The  relative  levels  of  transmission  loss  for  each  of  the 
four  environments  considered  here  remain  consistent  at  each 
receiver  depth.  In  general,  the  effect  of  shifting  the 
source  half  a  convergence  zone  from  0  to  30  km  is  negligable. 
In  fact,  the  transmission  loss  curves  for  these  two  cases 
often  resemble  closely  that  for  a  range-independent  environ¬ 
ment.  Consequently,  for  the  remainder  of  this  analysis, 
only  the  two  most  extreme  cases  will  be  considered,  viz., 
the  range-dependent  environment  with  the  source  located  at 
the  eddy  center  and  the  totally  range-independent  track. 
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7.4  Cold  Eddy  300  Hz 


An  examination  of  the  differences  in  transmission 
loss  for  these  same  source-receiver  geometries  at  300  Hz  reveals 
some  interesting  findings.  For  the  receiver  located  at  95  m 
(Figure  24),  the  range-dependent  environment  with  the  source 
at  the  eddy  center  shows  2  -  3  dB  more  loss  than  the  range- 
independent  environment.  This  is  the  reverse  of  what  was 
seen  earlier  for  the  cold  eddy  case  at  20  Hz.  This  is 
somewhat  of  an  anomalous  observation,  as  this  reversal  of 
the  relationship  exists  only  for  this  case  and  is  seen  in  no 
other  instance. 

Moving  to  a  receiver  depth  of  450  m,  (Figure  25) 
it  may  be  seen  that  this  relationship  has  reversed  itself, 
and,  once  again,  the  range-dependent  environment  with  the 
source  at  the  eddy  center  shows  consistently  less  loss  than 
the  range-independent  environment.  Figures  26  and  27, 
showing  transmission  loss  for  receivers  located  at  700 
and  900  m,  respectively,  echo  these  findings.  In  general, 
the  difference  between  the  two  environments  is  less  dramatic 
at  300  Hz  than  at  20  Hz. 


7.5  Warm  Eddy  300  Hz. 

Propagation  through  a  warm-core  eddy  at  300 
Hz  shows  similarly  consistent  findings.  For  a  receiver 
located  at  95  m  (Figure  28),  the  range-dependent  environment 
with  the  source  at  the  eddy  center  generally  shows  more  loss 
than  does  the  range-independent  environment,  with  the 
exception  of  the  first  50  nm  of  the  track.  This  difference 
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is  understandable  when  we  look  at  Figure  29,  which  shows  the 
unaveraged  transmission  loss  plots  for  the  same  environment. 
Transmission  loss  for  the  source  at  the  eddy  center  shows 
the  presence  of  a  surface  duct  in  the  first  90  km  of  the 
track.  Energy  remains  trapped  in  this  duct:  thus,  the 
observed  transmission  loss  for  this  range-dependent  environ¬ 
ment  is  actually  lower  than  the  corresponding  range-indepen¬ 
dent  environment.  Figure  30,  which  shows  transmission  loss 
for  a  receiver  at  450  m,  indicates  that  transmission  loss  is 
consistently  higher  for  the  range-dependent  environment  with 
the  source  at  the  eddy  center.  It  is  clear  that  the  surface 
duct  does  not  extend  to  this  depth,  as  transmission  loss  is 
lower  for  the  source  at  the  eddy  center,  even  in  the  first  50 
nm  of  the  track.  Figures  31  and  32,  for  receiver  depths  at 
700  and  900  m,  echo  these  findings. 

7.6  Cold  Eddy  1000  Hz. 

Figure  33  shows  transmission  loss  at  1  kHz  through 
the  cold  eddy  for  a  receiver  depth  of  95  m  and  the  source  at 
the  eddy  center.  In  the  first  185  km  of  the  track,  the 
range-dependent  environment  shows  less  loss  than  does  the 
range-independent  environment.  This  trend  reverses  itself 
past  185  km,  with  the  range-dependent  environment  showing 
slightly  more  loss  by  the  end  of  the  track.  The  markedly 
higher  level  of  transmission  loss  at  1  kHz  compared  to  300 
Hz  is  due  primarily  to  volume  attenuation,  which  plays  a 
significant  role  at  1  kHz.  Over  the  length  of  a  370  km 
track,  volume  attenuation,  computed  by  Thorp's  Equation,  is 
approximately  25  dB.  Volume  attenuation  plays  a  smaller 
role  at  300  Hz,  amounting  to  only  5  dB  attenuation  over  a 
400  nm  track. 
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7.7  Frequency  Effects 


The  effect  of  variations  in  source  frequency  on 
transmission  loss  is  apparent  from  Figure  .?4,  which  shows 
transmission  loss  in  the  range-dependent  environment  for  50. 
300,  and  1000  Hz.  The  constant  offset  in  transmission  loss 
between  the  20  and  300  Hz  cases  is  primarily  attributable  to 
surface  image  interference.  At  low  frequencies,  for  a 
shallow  source  and/or  receiver,  several  paths  combine 
destructively  to  yield  higher  levels  of  transmission  loss 
than  would  be  seen  for  a  deeper  source  and  receiver. 

The  effect  of  volume  attenuation  on  the  1000  Hz 
case  is  demonstrated  by  the  dramatically  higher  levels  of 
transmission  loss  for  that  case.  Recalling  that  volume 
attenuation  at  1  kHz,  over  200  nm,  is  approximately  25  dR. 
it  is  evident  that  were  it  not  for  volume  attenuation,  the 
300  and  1000  Hz  transmission  loss  curves  would  be  quite 
similar.  Thus,  the  largest  differences  attributable  to 
frequency  dependence  are  seen  between  20  and  300  Hz.  Above 
300  Hz,  there  is  relatively  little  additional  impact  on 
transmission  loss  resulting  from  the  eddy. 


fe-Dependent  Envlronne 
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Section  fi 

COMPARATIVE  ACOUSTIC  STATISTICS 


To  quantify  more  precisely  the  differences  between 
the  range-independent  and  range-dependent  environments  which 
were  seen  in  the  previous  transmission  loss  estimates,  we 
have  constructed  a  measure  of  statistical  difference  for 
comparative  purposes.  This  measure  is  actually  a  relative 
intensity  difference,  which  is  found  by  determining  the 
difference  in  intensity  between  any  two  environments  rela¬ 
tive  to  the  intensity  in  a  reference  environment.  The 
relative  intensity  difference  may  be  computed  on  a  point-by¬ 
point  basis  from  the  intensities  of  the  acoustic  field 
within  a  fixed  range-depth  window.  Figure  35  shows  a 
comparison  of  the  range-dependent  environment  with  the 
source  at  the  eddy  center  relative  to  the  range-independent 
environment  at  20  Hz.  The  range  window  for  this  comparison 
is  35  -  55  km.  The  depth  window  is  from  the  surface  to 
150  m.  The  histogram  shows  the  relative  frequency  of 
occurrence  within  this  window  of  the  statistic  (II  -  I?)/I2. 
where  II  is  the  intensity  in  the  range-dependent  environ¬ 
ment  and  12  is  the  reference  Intensity  that  of  the  range-in¬ 
dependent  environment.  The  histogram  reveals  that  the 
ma.iority  of  the  observations  are  at  the  positive  end  of  the 
scale,  indicating  that  the  range-dependent  environment  with 
the  source  at  the  eddy  center,  II,  has  a  higher  intensity 
over  this  range-depth  interval  than  does  the  range-indepen¬ 
dent  environment,  12.  This  is  reflected  in  the  mean  intensi¬ 
ties,  expressed  in  dB  units,  for  each  environment.  The  mean 
intensity  for  the  range-dependent  environment  is  91.3  dP, 
whereas  the  mean  intensity  for  the  range-independent  track 
in  the  same  window  is  94.8  dP,  3.5  dB  less. 
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The  nature  of  this  relationship  changes  with 
horizontal  range,  as  is  seen  in  Figure  36,  which  shows  a 
similar  comparison  for  the  range  window  between  55  and  75 
km.  In  this  instance,  the  distribution  is  positively  skewed, 
indicating  that  the  range-independent  environment  has 
more  energy  in  it,  for  this  range-depth  window,  than  does 
the  range-dependent  environment.  This  is  also  evident  in 
the  mean  intensities  for  these  two  environments,  which  show 
that  the  loss  associated  with  the  range-dependent  environ¬ 
ment  is  higher  in  this  instance  than  the  loss  associated 
with  the  range-independent  environment.  In  a  similar 
manner.  Figures  37  and  38  examine  the  distribution  of  energy 
in  this  fixed  depth  window  for  the  next  two  20  km  range 
windows. 

Figure  39  shows  the  relative  intensity  difference 
at  20  Hz  for  a  warm  eddy  (II)  versus  a  range-independent 
environment  (12).  The  distribution  in  this  instance  is 
positively  skewed,  indicating  that  the  range-dependent 
environment  with  the  source  at  the  eddy  center  shows  higher 
levels  of  loss  than  does  the  range  independent  environment. 
This  is  consistent  with  the  transmission  loss  curves  exam¬ 
ined  earlier.  Further,  as  both  II  and  12  are  intensities 
greater  than  or  equal  to  0,  the  smallest  value  of  the 
relative  Intensity  difference  (II  -  12) / 12  is  -1.  It  is 
evident  from  the  histogram  and  from  the  mean  value  of  this 
statistic  that  most  of  the  points  approach  this  lower  limit 
of  -1 . 

Figure  40  shows  relative  intensity  differences  for 
the  next  20  km  section  of  track,  55  -  75  km.  In  this 
instance,  the  distribution  is  actually  bimodal,  although  the 
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♦  Figure  37 


KtLRTIVE  INTENSITY  DIFFERENCES  FGR 

COLD  EDDY  (20HZ} —CENTER  (II)  VS  RfiNGE-1  112) 


DEPTH—  0.00  TO  SOO.OO  FT 

RaNGE—  ^0.00  TO  SO.  DO  NMI 

N=  119- 


MERNdl)  = 

3.483  X 

10-}}  = 

104.580 

□3 

SDEV(Il)  = 

6.2IJ5  X 

10-  = 

102.044 

03 

MEaN(12)  = 

0. 02S  X 

10-  = 

126.040 

03 

S0EV(I2)  = 

0.040  X 

10-"  = 

124.030 

DB 

MERN  (- ]  =  623.  m4 
SDEV  -  874.317 


RELRTIVE  FREOUENCY  (PER  CENT) 

0  20  40  BO  80  100 


RELRTIVE  FREOUENCY  (PER  CENT) 

8-5 


11-12 


Fi^^ire-SS 


RELRTIVE  INTENSITT  OIFFERFN'CES  FflR 

COLD  EDDY  (20HZ) —CENTER  (ID  VS  RRNGE-I  (J2) 


DEPTH—  O.QO  TO  500.00  FT 

RANGE—  50.00  TO  60.00  NM! 

N=  189 


NERNDD  = 

28.717 

X 

10' 

SDEVdll  = 

20.514 

X 

10' 

MEAN (12)  = 

15.983 

X 

10' 

SDEV (12)  = 

25.693 

X 

10' 

MEAN  [- 


SDEV  (- 


llzlZ) 

12  J 

11-12-1 
12  J 


107.717 


174.845 


95.419 
95.858 
97. 700 
95.902 


RELATIVE  FREQUENCY  (PER  CENTl 
Q  20  40  BO  80 


RELATIVE  FREQUENCY  (PER  CENT) 


Fijiure  39 


RELRTIVE  INTENSITY  DIFFERENCES  FOR 

WARM  EDDY  (20  HZJ — RS=CENTER  (ID  fiND  RfiNG£-J  (12) 


DEPTH—  0.00 
RANGE—  20.03 
N-  133 

MERNlIl)  = 
SDEVdl)  = 
HERN (12)  - 
SDEV(I2)  = 

HERN  j  : 


TO  500.00 
TO  30.00 


0.523 

X 

10-  : 

112.313 

Or; 

0.852 

X 

10-  = 

110.G93 

JB 

31.531 

X 

10-  = 

34.391 

OB 

45.3SG 

X 

10-"  = 

93.336 

OB 

SDEV  ( )  =  O'-  342 


<-1.000 


RELATIVE  FREQUENCY 

20  40  r 


(PER  CENT) 

r  80 


>  3- 000 


20  40  60  80 

REI  flTfVF  FRFQflFMrY  fPFP  TFWT) 


11-12 


Figure  40 
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ma.iority  of  observations  indicate  that  the  ranfte-depenflent 
environment  has  more  energy  in  this  ranee-depth  window  than 
does  the  ranee-independent  environment.  Examination  of  the 
mean  intensity  values  in  dB  units  for  the  two  environments, 
however,  shows  that  the  mean  intensity  in  this  window  is 
actually  higher  for  the  range- i ndependen t  environment  than 
it  is  for  the  range-dependent  environment.  This  Illustrates 
the  potentially  misleading  conclusions  that  may  be  drawn 
from  a  simple  comparison  of  mean  intensity  over  an  opera¬ 
tional  window  of  interest.  Rather,  the  use  of  the  histo¬ 
grams  contributes  significantly  different  information,  which 
may  be  important  from  a  tactical  perspective.  Similarly, 
Figures  41  and  42  depict  the  relative  intensity  differences 
for  the  next  40  km  of  the  track. 

Relative  intensity  differences  within  a  cold 
eddy  at  300  Hz  may  be  seen  in  Figure  43.  The  depth  window 
remains  the  same  as  that  for  the  20  Hz  cases,  from  0  to  l.'io 
m.  The  range  window  is  considerably  smaller  in  this 
instance,  here  a  2  km  window  from  to  57  km  alone  the 
track.  The  histogram  indicates  that  within  this  ranee-depth 
window,  the  range-independent  environment  contains  more 
energy  than  does  the  range-dependent  environment.  This  is 
also  reflected  in  the  mean  intensities  for  the  individual 
windows  and  in  the  mean  value  of  the  statistic  (II  -  I?)/!?. 

This  relationship  changes  somewhat  in  the  next 
range  window  considered,  between  74  and  7fi  km  (Figure 
44).  While  the  mean  intensity  for  each  environment  as  well 
as  the  mean  value  of  the  relative  intensity  difference 
indicates  that  the  range-independent  environment  has  more 
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energy  in  it  than  does  the  range  dependent  environment,  the 
nature  of  this  relationship  has  shifted  somewhat,  as  is  seen 
in  the  histogram. 
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Relative  intensity  differences  for  a  warm  eddy  at 
300  Hz  may  be  seen  in  Figure  45.  The  histogram  indicates 
that  a  significant  number  of  observations  within  this 
range-depth  window  show  that  the  range-independent  environ¬ 
ment  has  more  energy  associated  with  it  than  does  the 
range-dependent  environment.  The  mean  relative  intensity 
difference  within  this  range  depth  window,  however,  indi¬ 
cates  that  the  range  dependent  environment  has  more  energy 
associated  with  it  than  does  the  range  independent  environ¬ 
ment.  This  is  another  example  of  the  misleading  conclusions 
that  may  be  drawn  from  the  use  of  a  single  summary  statistic 
Figure  46  shows,  for  the  range  window  74  to  76  km,  a 
reversal  of  this  relationship.  In  this  instance,  the 
range-dependent  environment  shows  significantly  higher 
intensity  than  the  range  independent  environment. 
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Figure  45 
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Section  9 
SUMMARY 


In  summary,  the  eddy  model  developed  by  Henrick 
has  been  examined  in  detail.  This  model  examines  the 
perturbations  to  density,  temperature,  and  sound  velocity 
structures  which  result  from  the  presence  of  an  eddy. 
Statistically,  the  modeled  density,  temperature,  and 
sound  velocity  profiles  agree  fairly  well  with  the  back¬ 
ground  conditions.  The  largest  deviations  from  the  observed 
conditions  occur  in  the  upper  600  m.  This  is  due  primarily 
to  the  lack  of  surface  expression  in  the  model,  as  well  as 
to  its  inability  to  model  irregular  profiles.  This  poses 
a  particular  problem  from  an  acoustic  perspective  as,  for 
ASW  depths  of  interest,  significant  differences  between  the 
observed  and  the  modeled  conditions  will  result.  A  model 
upgrade  is  in  order,  to  permit  the  modeling  of  irregular 
profiles  and  eddy  surface  expressions.  Cubic  spline  tech¬ 
niques  have  been  suggested  as  one  means  for  modeling  the 
irregular  profiles. 

From  the  acoustic  analysis,  it  has  been  demonstra¬ 
ted  that  little  difference  exists  between  a  range-dependent 
environment  with  the  source  at  0  km  and  one  in  which  the 
source  is  located  at  30  km.  Shifting  the  source  by  half  a 
convergence  zone,  therefore,  exerts  relatively  little 
influence  on  transmission  loss.  The  most  significant 
differences  occurred  between  the  range-dependent  environment 
with  the  source  at  the  eddy  center  and  the  range-independent 
environment.  At  20  Hz,  within  a  cold  eddy,  the  range- 
dependent  environment  may  show  as  much  as  10  db  less  loss 
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than  the  range-independent  environment.  In  the  warm  eddy 
this  relationship  is  reversed,  with  the  range-dependent 
environment  showing  more  loss  than  the  range-indenendent 
environment.  At  300  Hz,  within  the  cold  eddy,  the  same 
general  trends  as  those  seen  at  20  Hz  are  observed,  with  the 
exception  of  95  m,  where  the  range-independent  environment 
shows  less  loss  than  the  range-dependent  environment.  For 
the  warm  eddy,  the  same  trends  are  observed  at  300  Hz 
that  were  observed  at  20  Hz.  In  general,  there  is  less 
overall  difference  between  the  two  environments  at  300  Hz 
than  at  20  Hz.  At  1000  Hz,  in  a  cold  eddy,  the  range- 
dependent  environment  shows  less  loss.  The  major  difference 
in  transmission  loss  between  300  Hz  and  1000  Hz  is  attribu¬ 
table  to  volume  attenuation. 

The  differences  observed  in  the  transmission 
loss  estimates  were  examined  in  more  detail  with  the  use 
of  comparative  statistics  such  as  a  relative  intensity 
difference.  These  measures  of  statistical  difference  depict 
the  changing  distribution  of  energy  as  a  function  of  range 
and  depth  along  a  given  track,  and  may  be  used  to  pinpoint 
more  precisely  areas  where  the  presence  of  an  eddy  may  have 
a  significant  impact  on  tactical  systems. 


Section  10 
CONCLUSIONS 


From  the  work  described  here,  three  principal 
conclusions  can  be  drawn: 

•  ft'ith  respect  to  ASW  applications,  eddies  have 
an  appreciable  influence  on  the  distribution 
of  energy  in  range  and  depth  only  when  the 
target  is  near  the  eddy  center.  Furthermore, 
the  influence  appears  greater  at  low  frequencies 
(20  Hz)  than  at  high  frequencies  (300  Hz)  and 
may  amount  to  10  dB,  on  average,  for  some 
geometries. 

•  The  Henrick  eddy  model  is  not  sufficiently 
flexible  in  its  present  form  to  fit  measured 
Pacific  data  with  enough  accuracy  to  be  used 
for  acoustic  studies.  However,  some  possible 
improvements  appear  to  be  straightforward. 

•  Some  progress  has  been  made  here  toward  devel¬ 
oping  an  objective,  useful  measure  of  the 
acoustic  influence  of  mesoscale  features. 
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